Theory of tunneling conductance of graphene NIS junctions 
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We calculate the tunneling conductance of a graphene normal metal-insulator-superconductor 
(NIS) junction with a barrier of thickness d and with an arbitrary voltage Vo applied across the bar- 
rier region. We demonstrate that the tunneling conductance of such a NIS junction is an oscillatory 
function of both d and Vo- We also show that the periodicity and amplitude of such oscillations 
deviate from their universal values in the thin barrier limit as obtained in earlier work [Phys. Rev. 
Lett. 97, 217001 (2006)] and become a function of the applied voltage Vo. Our results reproduces 
the earlier results on tunneling conductance of such junctions in the thin [Phys. Rev. Lett. 97, 
217001 (2006)] and zero [Phys. Rev. Lett. 97, 067007 (2006)] barrier limits as special limiting 
cases. We discuss experimental relevance of our results. 

PACS numbers: 74.45+c, 74.78.Na 
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I. INTRODUCTION 

Graphene, a two-dimensional single layer of graphite, 
has been recently fabricated by Novoselov et. al 1 . This 
has provided an unique opportunity for experimental ob- 
servation of electronic properties of graphene which has 
attracted theoretical attention for several decades 2 -. In 
graphene, the energy bands touch the Fermi energy at 
six discrete points at the edges of the hexagonal Bril- 
louin zone. Out of these six Fermi points, only two are 
inequivalent; they are commonly referred to as K and 
K' points 3 . The quasiparticle excitations about these 
K and K' points obey linear Dirac-like energy disper- 
sion. The presence of such Dirac-like quasiparticlcs is 
expected to lead to a number of unusual electronic prop- 
erties in graphene including relativistic quantum Hall ef- 
fect with unusual structure of Hall plateaus^. Recently, 
experimental observation of the unusual plateau struc- 
ture of the Hall conductivity has confirmed this theo- 
retical prediction 5 -. Further, as suggested in Ref. @, the 
presence of such quasiparticles in graphene provides us 
with an experimental test bed for Klein paradox. 7 

Another, less obvious but nevertheless interesting, con- 
sequence of the existence Dirac-like quasiparticles can be 
understood by studying tunneling conductance of a nor- 
mal metal-superconductor (NS) interface of graphene 8 . 
Graphene is not a natural superconductor. However, 
superconductivity can be induced in a graphene layer 
in the presence of a superconducting electrode near it 
via proximity effecl£&i£ or by possible intercalation with 
dopant molecules 11 . It has been recently predicted 8 that 
a graphene NS junction, due to the Dirac-like energy 
spectrum of its quasiparticles, can exhibit specular An- 
dreev reflection in contrast to the usual retro reflection 
observed in conventional NS junction al 13 . Such specu- 
lar Andreev reflection process leads to qualitatively dif- 
ferent tunneling conductance curves compared to con- 
ventional NS junctions 8 . The effect of the presence of 
a thin barrier region of thickness d — > created by ap- 
plying a large gate voltage Vq — > oo ( such that Vod is 



finite) between the normal and the superconducting re- 
gion has also been studied in Ref. It has been shown 
that in this thin barrier limit, in contrast to all normal 
metal-insulator-superconductor (NIS) junctions studied 
so far, the tunneling conductance of a graphene NIS junc- 
tion is an oscillatory function of the dimcnsionless barrier 
strength x — Vod/(tiVF), where vf denotes the Fermi ve- 
locity of graphene, with periodicity n. Further, it has 
also been demonstrated that the tunneling conductance 
reaches its maxima of 2Gq for \ — i n + 1/2)tt, where n 
is an integer. The latter result was also interpreted in 
terms of transmission resonance property of the Dirac- 
Bogoliubov quasiparticles 5 -. However, no such studies 
have been undertaken for NIS junctions with barriers of 
arbitrary thickness d and barrier potential Vo . 

In this work, we extend the analysis of Ref. and cal- 
culate the tunneling conductance of a graphene NIS junc- 
tion with a barrier of thickness d and with an arbitrary 
voltage Vo applied across the barrier region. The main 
results of our work are the following. First, we show that 
the oscillatory behavior of the tunneling conductance is 
not a property of the thin barrier limit, but persists for 
arbitrary barrier width d and applied gate voltage Vb, 
as long as d <C £, where £ is the coherence length of 
the superconductor. Second, we demonstrate that the 
periodicity and amplitude of these oscillations deviate 
from their values in the thin barrier limit and becomes 
a function of the applied voltage Vq. We point out that 
the barriers which can be realistically achieved in current 
experimental setups 5 do not necessarily fall in the thin 
barrier regime which necessitates a detailed study of arbi- 
trary barriers as undertaken here. Finally, we show that 
our analysis correctly reproduces the tunneling conduc- 
tance for both zero barrier— and thin barrier 14 as limiting 
cases. 

The organization of the rest of the paper is as follows. 
In Sec. HH wc develop the theory of tunneling conduc- 
tance for a barrier of thickness d <C £ and with a voltage 
Vb applied across the barrier region and demonstrate that 
they correctly reproduce the results of Refs. [3 andU as 
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limiting cases. The results obtained from this theory is 
discussed in Sec.lIIIl Finally, in Sec. HVl we discuss possi- 
ble experiments that can be performed to test our theory. 



II. CALCULATION OF TUNNELING 
CONDUCTANCE 

Let us consider a NIS junction in a graphene sheet oc- 
cupying the xy plane with the normal region occupying 
x < —d for all y as shown schematically in Fig. [TJ The 
region I, modeled by a barrier potential V , extends from 
x = — d to x = while the superconducting region oc- 
cupies x > 0. Such a local barrier can be implemented 
by either using the electric field effect or local chemical 
doping^. The region x > is to be kept close to an 
superconducting electrode so that superconductivity is 
induced in this region via proximity effect^. In the rest 
of this work, we shall assume that the barrier region has 
sharp edges on both sides. This condition requires that 
d ^ X = 2n/kp, where kp and A are Fermi wave- vector 
and wavelength for graphene, and can be realistically cre- 
ated in experiments 6 . The NIS junction can then be 
described by the Dirac-Bogoliubov-de Gennes (DBdG) 
equations^ 



H a -E F + U(r) 
A*(r) 



E, 



A(r) 
U(r) 



7~L a 



lp a = Elp a . 



(1) 



Here, i\i a = (ip A a,ipBa,ip*Aa'-*P*Ba) are the 4 component 
wavefunctions for the electron and hole spinors, the in- 
dex a denote K or K' for electron/holes near K and K' 
points, a takes values K'(K) for a = K(K'), Ep denote 
the Fermi energy which can be made non-zero either by 
doping or by applying a potential to the graphene sheet, 
A and B denote the two inequivalent sites in the hexago- 
nal lattice of graphene, and the Hamiltonian TL a is given 
by 



H a = -ihvp (a x d x + sgn(a)a v d y ) . 



(2) 



In Eq.[21 vp denotes the Fermi velocity of the quasiparti- 
cles in graphene and sgn(et) takes values ± for a — K(K'). 
The pair-potential A(r) in Eq. Q] connects the electron 
and the hole spinors of opposite Dirac points. We have 
modeled the pair-potential as 



A(r) = Aq exp(i 



(3) 



where Ao and 4> are the amplitude and the phase of the 
induced superconducting order parameter respectively 
and 6(x) denotes the Heaviside step function. 

The potential t/(r) gives the relative shift of Fermi en- 
ergies in normal, insulating and superconducting regions 
of graphene and can be modeled as 



U(v) 



-U 6(x) + V 0(-x)6(x + d). 



(4) 
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FIG. 1: (Color online) A schematic sketch of a graphene NIS 
junction. The dashed region sees a variable gate (shown as red 
filled region) voltage Vb which creates the barrier. Additional 
gate voltage Uo, which may be applied on the superconducting 
side, and the current source is not shown to avoid clutter. 



The gate potential Uo can be used to tune the Fermi sur- 
face mismatch between the normal and the superconduct- 
ing regions. Notice that the mean-field conditions for su- 
perconductivity are satisfied as long as Ao <§C (Uq + Ef); 
thus, in principle, for large Uq one can have regimes where 
A > Ej&. 

Eq. [T] can be solved in a straightforward manner to 
yield the wavefunction ip in the normal, insulating and 
the superconducting regions. In the normal region, for 
electron and holes traveling the ±x direction with a 
transverse momentum k y — q and energy e, the (un- 
renormalized) wavefunctions are given by 

ipff = (l,±e ±iQ ,0,0)exp[i(±fc„ 



sin(a) 



(l,±e 

(0,0,1,T' 
hvFq 
e + E F ' 



sin(a 



exp [i (±k' n x - 

hvpq 
e-Ep' 



qy)} i 

h QV)] , 



(5) 



where the wave- vector k n (k' n 
wavefunctions are given by 



for the electron (hole) 



k n (k n ) — 



'e+{-)E F 



q 2 , 



(6) 



and a(a') is the angle of incidence of the electron (hole). 
In the barrier region, one can similarly obtain 

= (l,±e ±i9 ,0,0)exp[i(±k b x + qy)}, 

$t = (0,0,l,Te ±ie ')exp[i(±^ + ^)], (7) 

for electron and holes moving along ±x. Here the angle of 
incidence of the electron(hole) 0{O') and the wavevector 
kb{k' b ) are given by is 

sin [9{9')] = hvpq/[e + (-)(Ep-V )], 



'e + (-)(E F -V ) 



(8) 
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Note that Eq. [7] ceases to be the solution of the Dirac 
equation (Eq. [T|) when E F — V and e = 0. For these 
parameter values, Eq. [1] in the barrier region becomes 
HaipB = which do not have purely oscillatory solutions. 
For the rest of this work, we shall restrict ourselves to the 
regime Vq > E F . 

In the superconducting region, the BdG quasiparticles 
are mixtures of electron and holes. Consequently, the 
wavefunctions of the BdG quasiparticles moving along 
±x with transverse momenta q and energy e, for (Uq + 
E F ) 3> A , e, has the form 



' eT i/3 )Te ±i( 7 -/3) )e -i0 )Te i(± 7 - 



x exp [i (±k s x + qy) — kx] 
sin (7) = hv F q/(E F + U ), 



(9) 



where 7 is the angle of incidence for the quasiparticles. 
Here the wavevector k s and the localization length 
can be expressed as a function of the energy e and the 
transverse momenta q as 



[{U + E F ) /hv F f - 

(hv F ) 2 k s 
[{U + E F )A sm(P)y 



(10) 



where (3 is given by 



/3 = cos-^e/Ao) if|e| < A , 
= -icosh^ 1 (e/A ) if|e|>A . 



(11) 



Note that for |e| > A , k becomes imaginary and the 
quasiparticles can propagate in the bulk of the supercon- 
ductor. 

Next we note that for the Andreev process to take 
place, the angles 9, 9' and a' must all be less than 90°. 
This sets the limit of maximum angle of incidence a. 
Using Eqns. \5\ and one finds that the critical angle of 
incidence is 



a c = ai 1) 9(V ~2E F ) + a^9(2E F -V ) 



(2)z 



,(2) = 



arcsin [|e — E F \ / (e + E F )] , 
wccsin[\e-\E F -V \\/(e + E F )}. (12) 



Note that in the thin or zero barrier limits treated in 
Refs. [H] andd, a c = for all parameter regimes. 

Let us now consider a electron-like quasiparticle inci- 
dent on the barrier from the normal side with an energy 
e and transverse momentum q. The basic process of or- 
dinary and Andreev reflection that can take place at the 
interface is schematically sketched in Fig. [2l As noted in 
Ref.[H, in contrast to conventional NIS junction, graphene 
junctions allow for both retro and specular Andreev re- 
flections. The former dominates when e, Ao E F so 
that a — —a' (Eq. [5]) while that latter prevails when 
E F <C e, A with a = a'. Note that in Fig. [21 we have 
chosen ta to denote a retro Andreev reflection for illus- 
tration purposes. In practice, ta includes both retro and 





FIG. 2: (Color online) A schematic sketch of normal reflection 
(r), Andreev reflection (ta) and transmission processes (t and 
t') at a graphene NIS junction. Note that in this schematic 
picture, we have chosen ta to denote a retro Andreev reflec- 
tion for illustration purpose. In practice, as discussed in the 
text, ta takes into account possibilities of both retro and spec- 
ular Andreev reflections. The electron and hole wavefunctions 
inside the barrier region is not sketched to avoid clutter. 



specular Andreev reflections. In what follows, we shall 
denote the total probability amplitude of Andreev reflec- 
tion as ta which takes into account possibilities of both 
retro and specular Andreev reflections. 

The wave functions in the normal, insulating and su- 
perconducting regions, taking into account both Andreev 
and normal reflection processes, can then be written as 13 



N 



■ nip 



h- 
N ' 



(13) 



where r and ta are the amplitudes of normal and An- 
dreev reflections respectively, t and t' are the amplitudes 
of electron-like and hole-like quasiparticles in the super- 
conducting region and p, q, m and n are the amplitudes 
of electron and holes in the barrier. These wavefunctions 
must satisfy the appropriate boundary conditions: 



^N\x=-d = *B\x=-d, *BU=0 = *SU=0- (14) 



Notice that these boundary conditions, in contrast their 
counterparts in standard NIS interfaces, do not impose 
any constraint on derivative of the wavefunctions at the 
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boundary. These boundary conditions yield 
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id 


= pe~ lkbd + qe ik " d , 
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—ia^ikj 


id 


= -pe ii - e ~ kbd "> + qe~ i{e ~ kbd \ 






t d 










= me-< d + ne< d , 


r A e 


-iot' g^fej 


t d 


= -me^'- k '» d )+ne- i( - e '- k '"\ 




p + 


<1 




-pe w 


+ qe~ 


iO 


= _ te i(7-« +t ' e -i(7-/3) s 




m + 


n 




-me 10 ' 


+ ne~ l 


8' 


= -te^-V +t , e -*(T+^. (15) 



Using the boundary conditions fEq. II 5[) . one can now 
solve for the coefficients r, r A , t and t' in Eq. [T31 After 
some straightforward but cumbersome algebra, we find 
that 



-2ifc„d 



(16) 



Af = [e ia cos(k b d + 6)-i sm{k b d)] 

-p[cos(k b d -9)-i e ia sin(fc b d)] , (17) 
V = [e~ ia cos(M + 6) + i sin(fc 6 d)] 

+p [cos(M - 0) + ie- ia sin(M)] , (18) 

= cos(g)[Te-^ + W\ V COs( - hd ^ sin (M] 

+re lfc " ti [cos(fc f) d -9) + ie- ia sin(fc 6 <f)]) , (19) 
t = Tt', (20) 
<(r + l)e< d cos(6»')e-^ 



TA 



cos(k' b d - 9') - ie~ ia ' sm(k' b d) ' 



(21) 



where the parameters T and p can be expressed in terms 
of 7, (3, 6, 6', a, and a' (Eqs. El 1 H andrjl]) as 



r 
v 



_p e »(7-^) 4. e -i(7-/3) 



-27 



e" ia ' cos(fc£d + 6') - i sm(k' b d) 
cos(k' b d - 9') - ie~ ia ' sm(k' b d) ' 



(22) 
(23) 
(24) 



The tunneling conductance of the NIS junction can now 
be expressed in terms of r and ta byi^ 



GjeV) = 
Go (eV) J Q 



a,, ,2 cos(a')\ 
1 — \r\ + \ta\ r ^)cos(a)aa, 



cos(a) 



(25) 



where Go = 4e 2 iV(eU)//i is the ballistic conductance 
of metallic graphene, eV denotes the bias voltage, and 
N(e) = (Ep + e)w/(irhvp) denotes the number of avail- 
able channels for a graphene sample of width w. For 
eV <C Ep, Go is a constant. Eq. [23 can be evaluated 
numerically to yield the tunneling conductance of the 
NIS junction for arbitrary parameter values. We note 



at the outset, that G = when a c — 0. This occurs 
in two situations. First, when eV — Ep and Vo > 2Ep 

so that a c — ai^ vanishes. For this situation to arise, 
Ep + Uo > A > Ep which means that Uo has to be finite. 
Second, a c — a y c ' — when eV = and Ep = Vo, so 
that the zero-bias conductance vanishes when the barrier 
potential matches the Fermi energy of the normal side^. 

We now make contact with the results of Ref. [lj in the 
thin barrier limit. We note that since there are no condi- 
tion on the derivatives of wavefunctions in graphene NIS 
junctions, the standard delta function potential approx- 
imation for thin barrier— can not be taken the outset, 
but has to be taken at the end of the calculation. This 
limit is defined as d/\ — * and Vo/Ep — > 00 such that 
the dimensionless barrier strength 



X = Vod/hvp = 2ir 



Vo_ 

E F 



(26) 



remains finite. In this limit, as can be seen from Eqs. [5] 
|8]and[9l 9, 9', k n d, k' n d — » and k b d, k' b d ->xso that the 
expressions for T, p and r\ (Eq. l24jl 



r 



tb 



e -*7 _ ^tb 
e i7 + jytb 



<1 



,tb 



e ia ' cos(x) — i sin(x) 
cos(x) — ie~ la ' sin(x) 



„tb 



e -i{l-0) _ ptb e i( 7 -/3) 
ptbg— i/3 _|_ gi/3 



(27) 



where the superscript "tb" denotes thin barrier. Using 
the above-mentioned relations, we also obtain 



.,1b 



'tb 



.tb 



7 A 



cos(x) (e ia - p th ) - i sin(x) (l - p th e la ) 
cos(x) (e~ M + p tb ) + i sin(x) (1 + p tb e~ ia ) ' 



ix) (1 



^tb'i 



3in(x) (e 



Tt' th , 

t nb ( r + X ) e -i<t> 

cos(x) — ie~ la ' sin(x) 



(28) 



Eq s. l27l and [28l are precisely the result obtained in Ref. 
Il4l for the tunneling conductance of a thin graphene NIS 
junction. The result obtained in Ref. [8j can be now easily 
obtained from Eqs. [27] and [28] by substituting % = in 
these equations, as also noted in Ref. 



III. RESULTS 

A. Qualitative Discussions 

In this section, we shall analyze the formulae for tun- 
neling conductance obtained in Sec. [Tl] First we aim to 
obtain a qualitative understanding of the behavior of the 
tunneling conductance for finite barrier strength. To this 
end, we note from Eq. [25] that the maxima of the tunnel- 
ing conductance must occur where |r| 2 is minimum. In 
fact, if |r| 2 = for all transverse momenta, the tunneling 
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FIG. 3: (Color online) Plot of zero-bias tunneling conductance 
for Uo = and Ao = 0.01-Ef as a function of gate voltage Vb 
and barrier thickness d. Note that the oscillatory behavior of 
the tunneling conductance persists for the entire range of Vb 
and d. 



conductance reaches its value 2Go- Therefore we shall 
first try to analyze the expression of r (Eq. I16|) for sub- 
gap voltages and when the Fermi surfaces of the normal 
and superconducting sides are aligned with each other 
(Uo = 0). In this case, we need A <C Ep. So for subgap 
tunneling conductance, we have e < A <C Ep. In this 
limit, a ~ -a' ~ 7 (Eqs. [5]and[pJ, k b ~ k' b , and 9 ~ -0' 
(Eq. [5]). Using these, one can write 



50 

V o /E p 



75 



100 



V 



P 



e ia cos(k b d — 9) — i sin(fc&d) 
cos(k b d + 9) — ie ta sm(k b d) ' 
rj cos(a — 0) + i sin(/3) 
cos(a + j3) + ir)sm(P) 



(29) 
(30) 



Substituting Eq. [3D] in the expression of TV, we find 
that the numerator of the reflection amplitude r becomes 
(Eqs. Dl and EH) 



A/" = t— 



e' 
D 



4sin(a) sin(/3) cos(k b d — 9) 



x 



i cos(a) sm(kbd) 

+ (cos(k b d -0)+ cos(k b d + 6))/2 
+2 [cos(k b d + 8) - cos(k b d - 0)} 
x cos(a — (3) {cos(a) + [cos(k b d — 9) 

+ cos(k b d + 9)} /2} + sm(k B d) sin(/3) 



(31) 



Do = cos(k b d + 9) cos(a + 0) + sin(k b d) sin(/3) 

+ie la [cos(k b d — 9) sin(/3) — sin(k b d) cos(a + 0)] . 

(32) 

From the expression of TV (Eq.GH]), we n °te the follow- 
ing features. First, for normal incidence (a = 0) where 



FIG. 4: Plot of periodicity Xpcriod of oscillations of tunneling 
conductance as a function of applied gate voltage Vb for Uo = 
and Ao = 0.01-Ef- Note that the periodicity approaches 
7r as the voltage increases since the junction approaches the 
thin barrier limit. 



9 = 9' = 0, M and hence r (Eq. [THJ) vanishes. Thus 
the barrier is reflectionless for quasiparticles which inci- 
dent normally on the barrier for arbitrary barrier thick- 
ness d and strength of the applied voltage Vb- This is 
a manifestation of Klein paradox for Dirac-Bogoliubov 
quasiparticles^. However, this feature is not manifested 
in tunneling conductance G ( Eq. I25p which receives con- 
tribution from all angles of incidence. Second, apart from 
the above-mentioned cases, r never vanishes for all an- 
gles of incidence a and arbitrary eV < Ao unless 9 = 9'. 
Thus the subgap tunneling conductance is not expected 
to reach a maximum value of 2Go as long as the thin 
barrier limit is not satisfied. However, in practice, for 
barriers with Vb > 4Ep, the difference between 9 and 9' 
turns out to be small for all q < hp (< 0.25 for q < hp 
and eV = 0) so that the contribution to M (Eq. |3"T]) 
from the terms ~ (cos(k b d + 9) — cos(k b d — 9)) becomes 
negligible. Thus |r| 2 can become quite small for spe- 
cial values of Vb for all q < kp so that the maximum 
value of tunneling conductance can reach close to 2Go- 
Third, for large Vb, for which the contribution of terms 
~ (cos(k b d + 9) — cos(k b d — 9)) becomes negligible, Af 
and hence r becomes very small when the applied volt- 
age matches the gap edge i.e. sin(/3) = (Eq. l3"Tj) . Thus 
the tunneling conductance curves approaches close to its 
maximum value 2Gq and becomes independent of the 
gate voltage Vb at the gap edge eV = A for A <C Ep, 
as is also seen for conventional NIS junctions^. Fourth, 
in the thin barrier limit, (Vo/Ep — * 00 and d/X — * 0), 
9 — > and k b d — > x, so that the contribution of the 
terms ~ (cos(k b d + 9) — cos(k b d — 9)) in Eq . |3"T1 vanishes 
and one gets 
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FIG. 5: Plot of the amplitude [G max (eV = 0; V )-G min (eV = 
0; Vo)]/Go = (G m ax — G m in)/Go of zero-bias tunneling conduc- 
tance as a function of the applied gate voltage Vb for Uq = 
and Ao = 0ME F . Note that G reaches 2G for V > AE F 
where the amplitude become independent of the applied gate 
voltage as in the thin barrier limit and vanishes for Vo/E F = 1 
as discussed in the text. 



tb _ 2 sin(a) [sinQ + (3) - sin(x - (3)} 

K 

x [— cos(x) + i sin(x) cos(a)] , (33) 
£>o b = cos(x) cos(a + (3) + sin(x) sin(/3) + ie ia 

x [cos(x) sin(/3) - sin(x) cos(a + 0)] . (34) 



As noted in Ref. Q M tb and hence r th (Eq. |25J) vanishes 
at x = { n + 1/2)tt which yields the transmission reso- 
nance condition for NIS junctions in graphene. Fifth, 
as can seen from Eqs. HH] and HU both |r| 2 and \ta\ 2 
are periodic functions of Vq and d since both kf, and 
8 depend on Vq. Thus the oscillatory behavior of sub- 
gap tunneling conductance as a function of applied gate 
voltage Vo or barrier thickness d is a general feature 
of graphene NIS junctions with d <C £. However, un- 
like the thin barrier limit, for an arbitrary NIS junc- 
tion, k b d = xV( e f/Vo - l) 2 + hWpqt/Vg # x, and 
9^0. Thus the period of oscillations of |r| 2 and Ir^j 2 
will depend on Vq and should deviate from their univer- 
sal value 7T in the thin barrier limits. Finally, we note 
from Eqs. HH [25] and [33] that in the thin barrier limit 
(and therefore for large Vb), the amplitude of oscilla- 
tions of the zero-bias conductance for a fixed Vb, defined 
as [G max (eV = 0;Vq) - G min (eV = 0;Vb)]/G O) which 
depends on the difference of \r{\ = (n + l/2)7r)| 2 and 
\ r (x — nir)\ 2 becomes independent of \ or the applied 
gate voltage Vq. 



FIG. 6: (Color online) Plot of tunneling conductance as 
a function of the bias-voltage eV and gate voltage Vo for 
d = 0.4A and A = 0.01E F . Note that for large Vo, the tun- 
neling conductance at eV = Ao is close to 2Go and becomes 
independent of Vo (see text for discussion). 



B. Numerical Results 

The above-mentioned discussion is corroborated by nu- 
merical evaluation of the tunneling conductance as shown 
in Figs. [3] S][5] and [S] From Fig. [3] which plots zero-bias 
tunneling conductance G(eV — 0) as a function of Vo and 
d, we find that G(eV — 0) is an oscillatory function of 
both Vo and d and reaches close to its maximum value of 
2Go throughout the plotted range of Vb and d. Further, 
as seen from Fig. [4] the periodicity of these oscillations 
becomes a function of Vo- To measure the periodicity of 
these oscillations, the tunneling conductance is plotted 
for a fixed Vb as a function of d. The periodicity of the 
conductance <i per i d is noted down from these plots and 
Xperiod = Vodporiod/^F is computed. Fig. 0]clearly shows 
that Xperiod deviate significantly from their thin barrier 
value 7r for low enough Vb and diverges at Vq — > Ep 16 . 
Fig. [5] shows the amplitude of oscillations of zero-bias 
conductance as a function of Vb . We note that maximum 
of the zero-bias tunneling conductance G max (eV = 0) 
reaches close to 2G for V > V 0c ~ AE F . For V > V 0c , 
the amplitude becomes independent of the applied volt- 
age as in the thin barrier limit, as shown in Fig. [5] For 
Vo E F , a c = a { c ] -> 0, so that G(eV = 0) -> and 
hence the amplitude vanishes. Finally in Fig. [5] we plot 
the tunneling conductance G as a function of the applied 
bias-voltage eV and applied gate voltage Vb for d = 0.4A. 
We find that, as expected from Eq. [33] G reaches close 
to 2Go at the gap edge for all Vb > 6E F - Also, as in 
the thin barrier limit, the oscillation amplitudes for the 
subgap tunneling conductance is maximum at zero-bias 
and shrinks to zero at the gap edge eV = Aq, where the 
tunneling conductance become independent of the gate 
voltage. 

Next, we consider the case Uq ^ 0, so that Aq — E F <C 
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FIG. 7: (Color online) Plot of tunneling conductance as a 
function of the bias-voltage eV and the gate voltage Vb for 
d = 0.4A, A = 2E F and f/o = 25E F . As discussed in the 
text, the tunneling conductance is virtually independent of 
the applied gate voltage Vo due to the presence of a large Uo- 
Note that maximum angle of incidence for which Andreev 
reflection can take place vanishes at eV = Ep leading to 
vanishing of G at this bias voltage. 



(Ep + Uo)- In this regime, there is a large mismatch of 
Fermi surfaces on the normal and superconducting sides. 
Such a mismatch is well-known to act as an effective bar- 
rier for NIS junctions. Consequently, additional barrier 
created by the gate voltage becomes irrelevant, and we 
expect the tunneling conductance to become indepen- 
dent of the applied gate voltage V$. Also note that at 
eV = Ef, a c = (Eq. [T2|) . Hence there is no Andreev 
reflection and consequently Go vanishes for all values of 
the applied gate voltage for this bias voltage. Our re- 
sults in this limit, coincides with those of Ref. Finally 
in Fig. [5J we show the dependence of amplitude of oscil- 
lation of zero-bias tunneling conductance on Uo for the 
applied bias voltages Vo = 6Ep and Ao = O.QlEp. As 
expected, the oscillation amplitude with decreases mono- 
tonically with increasing Uq. We have verified that this 
feature is independent of the applied gate voltage Vq as 
long as V > V 0c . 



IV. EXPERIMENTS 

Superconductivity has recently been experimentally 
realized in grapheneii. In our proposed experiment, 
one needs to realize an NIS junction in graphenc. The 
local barrier can be fabricated using methods of Ref. 
5. The easiest experimentally achievable regime corre- 
sponds to Aq <C Ep with aligned Fermi surfaces for the 



FIG. 8: Plot of amplitude of oscillation (G max — G m i n )/Go 
of zero-bias tunneling conductance as a function of Uo/Ep 
for Vb = 6Ef and Ao = 0.01-Ef- The oscillation amplitudes 
always decay monotonically with increasing Uo independent 
ofVo. 



normal and superconducting regions. We suggest mea- 
surement of tunneling conductance curves at zero-bias 
(eV = 0) in this regime. Our prediction is that the 
zero-bias conductance will show an oscillatory behavior 
with the bias voltage. In graphene, typical Fermi en- 
ergy can be Ep < 40meV and the Fermi-wavelength is 
A > 100nm£&±!. Effective barrier strengths of < 80meV£ 
and barrier widths of d ~ 10 — 50 nm therefore specifies 
the range of experimentally feasible junctions^. Con- 
sequently for experimental junctions, the ratio Vo/Ep 
can be arbitrarily large within these parameter ranges 
by fixing Vo and lowering Ep. Experimentally, one can 
set 5 < Ep < 20meV so that the conditions Ao <C Ep 
Vo/Ep 3> 1 is easily satisfied for realistic Ao ~ 0.5meV 
and Vb = 200meV. This sets the approximate range 
Vo/Ep > 10 for the experiments. Note that since the 
period (amplitude) of oscillations increases (decreases) 
as Vo/Ep — * 1, it is preferable to have sufficiently large 
values of Vo/Ep for experimental detection of these os- 
cillations. 

To check the oscillatory behavior of the zero-bias tun- 
neling conductance, it would be necessary to change V 
in small steps SVo- For barriers of a fixed width, for 
example with values of d/X = 0.3, it will be enough to 
change Vb in steps of approximately 20 — 30meV, which 
should be experimentally feasible. We note that for the 
above-mentioned range of Vo/Ep, the experimental junc- 
tions shall not always be in the thin barrier limit. For 
example, as is clear from Fig. [4] the periodicity of oscil- 
lations Xpcriod of the zero-bias tunneling conductance of 
such junctions shall be a function of Vo and shall differ 
from 7T. This justifies our theoretical study of NIS junc- 
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tions in graphene which are away from the thin barrier 
limit. 

Apart from the above-mentioned experiments, it 
should also be possible to measure the tunneling conduc- 
tance as a function of the applied bias voltage eV/ Ao for 
different applied gate voltages Vq. Such measurements 
can be directly compared with Fig. [5J Finally, it might 
be also possible to create a relative bias Uq between the 
Fermi surfaces in the normal and superconducting side 
and compare the dependence of oscillation amplitudes of 
zero-bias tunneling conductance on Uq with the theoret- 
ical result shown in Fig. [7] 

In conclusion, we have presented a theory of tunnel- 
ing conductance of graphene NIS junctions with barriers 
of thickness d -C £ and arbitrary gate voltages Vq ap- 



plied across the barrier region. We have demonstrated 
that the oscillatory behavior of the tunneling conduc- 
tance, previously derived in Ref. for junctions with 
thin barriers, persists for all such junctions. However, 
the periodicity and amplitude of these oscillations devi- 
ate from their universal values in the thin barrier limit 
and become functions of the applied barrier voltage Vq. 
We have also shown that our work, which extends the 
earlier results of Ref. 0, correctly reproduce the earlier 
results for tunneling conductance obtained for thin^ and 
zerc£ barriers as limiting cases. We have discussed ex- 
perimental relevance of our results. 

KS and SB thank Graduate Associateship Program at 
Saha Institute which made this work possible. SB thanks 
T. Senthil and V.B. Shenoy for stimulating discussions. 
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